Systematic analysis of a spin-susceptibility representation of the pairing interaction in 

the 2D Hubbard model 
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A dynamic cluster quantum Monte Carlo algorithm is used to study a spin susceptibility represen- 
tation of the pairing interaction for the two-dimensional Hubbard model with an on-site Coulomb 
interaction equal to the bandwidth for various doping levels. We find that the pairing interaction 
is well approximated by |l7(r)^x(^ ~ -^0 with an efi'ective temperature and doping dependent 
coupling U{T) and the numerically calculated spin susceptibility x{K — K'). We show that at 
low temperatures, U may be accurately determined from a corresponding spin susceptibility based 
calculation of the single-particle self-energy. We conclude that the strength of the d-wave pair- 
ing interaction, characterized by the mean-field transition temperature, can be determined from a 
knowledge of the dressed spin susceptibility and the nodal quasiparticle spectral weight. This has 
important implications with respect to the questions of whether spin fluctuations are responsible 
for pairing in the high-Tc cuprates. 



I. INTRODUCTION 

Recent numerical calculations have shown that the 
dominant contribution to the d-wave pairing interaction 
in the 2D Hubbard model comes from the spin S = 1 
channel^ Motivated by this result, a simple spin sus- 
ceptibility representation of the pairing interaction was 
studied^. Results for a Hubbard on-site Coulomb inter- 
action equal to the bandwidth and a site filling (n) — 0.85 
have shown that the pairing interaction can be well ap- 
proximated by a simple RPA form^i^i^ii 



-U^{T)x{K-K'). 



(1) 



Here it was important that an effective temperature de- 
pendent coupling U{T) and the dressed spin suscepti- 
bility were used in Eq. ([T]) instead of the bare U and 
the perturbative RPA susceptibility. The coupling tJ{T) 
was determined by fitting the low frequency d-wave pro- 
jected irreducible particle-particle vertex calculated with 
a dynamic cluster approximation (DCA) quantum Monte 
Carlo (QMC) t echnique^t^ ' ^ ° i ^ ^ with the d-wave projec- 
tion of the form given by Eq. ([T]). Using this estimate 
of U{T) and the calculated dressed susceptibility and 
Green's function, it was shown that the eigenvalue and 
eigenfunction of the homogeneous Bethe-Salpeter equa- 
tion in the particle-particle channel are well represented 
by the corresponding quantities calculated with the ap- 
proximate interaction given by Eq. ([1]). 

Here we extend this study to explore this approxi- 
mation for other fillings and investigate other ways to 
determine the coupling strength U{T) that do not re- 
quire knowledge of the irreducible particle-particle ver- 
tex, since this is not experimentally accessible for the 
cuprates. We first discuss the DCA QMC technique used 
to calculate the relevant quantities and review the fitting 
procedure used to determine the temperature dependent 



coupling U{T) from the irreducible particle-particle ver- 
tex. We then examine how well Eq. ([T]) can describe 
the pairing interaction for various site fillings (n). We 
are particularly interested in the low doping regime. In 
this pseudogap opens in the density of states and 

the low-energy spin excitations^ ^i^'^i^^i^^i^^ and it is un- 
clear whether the form given by Eq. ([!]) is still a good 
representation of the pairing interaction. We then ex- 
plore two approximations that estimate U{T) from the 
single-particle spectrum. By assuming that the self- 
energy is determined by the same interaction [Eq. ([1])] 
as the particle-particle interaction, one can get a single- 
particle estimate for U{T). We first do this assuming 
that the exact single-particle spectral weight is known. 
Next we consider a scenario, where only limited infor- 
mation is available for the single-particle spectrum, such 
as the nodal spectral weight. In this case, one can use 
the interaction in Eq. ([T]) to self-consistently determine 
the self-energy and dressed Green's function in addition 
to the effective coupling U{T). This, of course, assumes 
that a reliable estimate of the non-interacting Green's 
function is available. 



We will study the quality of these approximations by 
comparing their respective estimates of the d-wave eigen- 
value with the "exact" result. Throughout this paper, 
"exact" will refer to the numerical results obtained using 
the DCA QMC technique. These DCA QMC calculations 
are carried out on a 4-site cluster. This means that phase 
fluctuations are suppressed and the temperature Tco at 
which the d-wave eigenvalue equals one corresponds to 
the mean- field transition temperature^^. This tempera- 
ture provides a natural measure of the strength of the 
d-wave pairing interaction. 
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II. DYNAMIC CLUSTER QUANTUM MONTE 
CARLO TECHNIQUE 

To calculate the single-particle self-energy, the pairing 
interaction and the spin susceptibility in the 2D Hub- 
bard model, we use a DCA QMC algorithuM:^ . The 
dynamic cluster approximation maps the original lattice 
model onto a periodic cluster of size Nc sites embed- 
ded in a self-consistent host. The essential assumption 
is that short-range quantities, such as the self energy 
and its functional derivatives (the irreducible vertex func- 
tions) are well represented as diagrams constructed from 
the coarse-grained Green's function. For the problem 
of interest, this is a reasonable assumption for systems 
where the correlations that mediate the pairing are short- 
ranged. To this end, the first Brillouin zone is divided 
into Nc cells, with each cell represented by its center 
wave-vector K surrounded by N/Nc lattice wavevectors 
labeled by k. The reduction of the TV-site lattice prob- 
lem to an effective Nc site cluster problem is achieved 
by coarse-graining the single-particle Green's function, 
i.e. averaging G(K + k) over the k within a cell which 
converges to a cluster Green's function Gc(K). Con- 
sequently, the compact Feynman diagrams constructed 
from G'c(K) collapse onto those of an effective cluster 
problem embedded in a host which accounts for the fluc- 
tuations arising from the hopping of electrons between 
the cluster and the rest of the system. The compact 
cluster quantities are then used to calculate the corre- 
sponding lattice quantities. 

The pairing interaction is given by the irreducible part 
of the particle-particle vertex 

TPP{K; K') = TPP{K, -K; K' , -K') (2) 

with K = (K,cj„). One can also use the DCA to calcu- 
late the spin susceptibihty x{Q,t^n)^'^- We then intro- 
duce a d-wave coupling strength^ 

l(.9(K)rgg,„(K, nT; K^ ^r)g(KO)KK' 

(3^(K))k ^ ' 

with the even frequency, even momentum part of the 
irreducible particle-particle vertex, 

rPP,„(K,7rr;K',7rr) = ^ ( rPP(K, ^T; K', ^T) 

-I- rPP(K,7rr,-K',7rr) 
-I- TPP{K,TrT,K',~TTT) 
+ rPP(K,7rr,-K',-7rr)) (4) 

and g(K) = (cos AT^; — cosAT^). The leading low temper- 
ature eigenvalue of the particle-particle Bethe-Salpeter 
equation is then calculated from 

„ (K, ~K- K', -K') xl^{K') MK') = 

\MK).{b) 



It is found to correspond to an eigenfunction with d-wave 
symmetryii2i. Here we have coarse-grained the Green's 
function legs, xf{K') = ^ Ek' Gt(K' + k') Gx(-K' - 
k'), according to the DCA assumption. Troughout this 
paper, we show results calculated on a 2x2 cluster for 
a near-neighbor hopping t — 1 and a Hubbard Coulomb 
interaction U — 8. 



III. SPIN SUSCEPTIBILITY 
REPRESENTATION 

InRef.i, we introduced an effective coupling strength 
U{T) by requiring that the d-wave coupling strength 
given by Eq. ^ be the same at a given temperature 
when TPP{K; K') is replaced by the approximate interac- 
tion given by Eq. ([1]) . 
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FIG. 1: (Color online) The coupling strength C/(T) versus 
temperature obtained from fitting the "exact" DCA QMC 
pairing interaction Eq. (|4]) for [/ = 8 and for diflterent values 
of the site filling (n). We work in units where t = \. 

In Fig. [1] we show the results for [/(T) for three dif- 
ferent fillings. For temperatures T < 1, [/(T) decreases 
with temperature for all fillings. One also sees that U de- 
creases with increasing doping for all temperatures. This 
result is consistent with earlier quantum Monte Carlo 
calculations that found that the electron-spin fluctua- 
tion vertex decreased with decreasing temperature and 
increasing doping^- . 

Using these estimates of C/(T) one can then explore 
how weh \U'^x{K - K') represents Tpp{K] K') by com- 
paring the d-wave eigenvalues. The curves with solid 
symbols in Fig. [5] show the d-wave eigenvalue versus T 
obtained from Eq. ([5]) with the "exact" DCA QMC in- 
teraction yPP . The curves with open symbols show the 
d-wave eigenvalue obtained from Eq. ([5|) when F^^ is re- 
placed by |c7^x(Ar — K'). In these calculations we have 
used DCA QMC results for x{K - K') as weU as the 
single-particle propagator G{k) that appears in Eq. 
With tj{T) determined from fitting F^^, we find for all 
fillings that the temperature dependence and size of the 
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FIG. 2: (Color online) The d-wave eigenvalue versus tem- 
perature obtained from the RPA form [Eq. ((T|] (open sym- 
bols) and from the "exact" DCA QMC interaction (solid sym- 
bols) for different values of the site filling {n). 



d-wave eigenvalue are reasonably accounted for by the 
simple form of the interaction given in Eq. ([T]). Since 
the DCA QMC calculations have been carried out us- 
ing an iVc = 4 site cluster, the phase fluctuations are 
suppressed and the temperature Tco at which the d-wave 
eigenvalue equals one represents a mean-field transition 
temperature. In Table H] we list the "exact" Tco value 
obtained when T^p is used to determine \d{T). The tem- 
perature TjIq^ at which = 1 when r^^* is replaced by 
3/2U^{T)x{K-K') with U{T) determined from Eq. ^ 
is also listed in Table HI One sees that this approximation 
over-estimates the mean-field transition temperature by 
of order 10-30%, depending upon the doping. Reasons 
for the disagreement at small doping could include the 
presence of the pseudogap in the spin excitations and the 
assumption of a frequency and K independent coupling 
strength U. 



(n) 



0.95 
0.90 
0.85 



Tco 



•.{1) 



T, 



(2) 



(3) 



0.080 0.100 (25%) 0.108 (35%) 0.105 (31%) 
0.074 0.087 (18%) 0.084 (14%) 0.081 (9%) 
0.067 0.074 (10%) 0.064 (4%) 0.058 (13%) 



TABLE I: The superconducting mean-field transition temper- 
ature obtained as the temperature where the d-wave eigen- 
value Atj = 1 for different values of the site filling (n) and 
different approximations. Tco'- "Exact" result obtained when 
the DCA QMC result for r"" is used; T^^^: result obtained 
when r''*' is replaced by the RPA form [Eq. H}] and U is de- 
termined from fitting the irreducible particle-particle vertex; 

(2) ^ 

T^Q ■ result otained when U is determined from fitting the 
nodel quasiparticle spectral weight using Eq. ®; rio^: result 
obtained when U and G{K) are determined self-consistently 
by fitting the nodal quasiparticle weight. The numbers in 
brackets denote the deviation from the "exact" result ex- 
pressed in percent. 



IV. SINGLE-PARTICLE FIT OF U{T) 

In the previous section, we determined the coupling 
strength U{T) by fitting the pairing interaction. Next we 
explore how well U{T) can be estimated from the single- 
particle spectrum by assuming that the self-energy ^{K) 
is determined by the same form, Eq. ([T]), i.e. given by 



;(72^Gc(i^-g)x(g) 



(6) 



with K = (K,ia;„). As was done for the pairing interac- 
tion and the Bethe-Salpeter equation [Eq. ([5|)], we first 
examine what happens if we use DCA QMC results for 
the susceptibility xiQ) ^-nd the sin gle- particle propaga- 
tor Gc{K). It was shown in Ref. ly that this simple 
representation of the self-energy provides a useful de- 
scription of the single-particle spectral weight A(k, w). 
Within this framework, we propose to estimate the cou- 
pling strength U{T) by requiring that the Matsubara 
quasiparticle weight 



Z(K,T) = 



1 - 



3mI](K,7rT) 



(7) 



for K = (7r/2,7r/2) calculated with the "exact" result for 
'2,{K) is the same at a given temperature for the approx- 
imate self-energy given by the form in Eq. 

For the 2x2 cluster, the DCA self-energy Y.{K) 
is calculated for the discrete set of momenta K = 
(0,0), (7r,0), (0,7r) and (7r,7r). To obtain the self- 
energy for k = (7r/2,7r/2) we interpolate T.{K) ac- 
cording to I](k, Wn) = e*'''^S](R, w„) where R 
are the distances in the cluster and S(R, a;„) is the 
Fourier-transform of S(K,a'„). This interpolation 
does not introduce fast Fourier components correspond- 
ing to length-scales larger than the cluster, and for 
the 2x2 cluster is guaranteed to preserve causality. 
For k — (7r/2,7r/2), one obtains I](7r/2, 7r/2, aj„) = 
\ [S(0, 0, c^„) + 2S(7r, 0, c^„) + I](7r, ^, w„)]- 

From a phenomenological point of view, it is interesting 
to see how well the d-wave eigenvalue calculated with the 
approximate pairing interaction in Eq. ([1]) and the cou- 
pling U determined from the self-energy reproduces the 
"exact" d-wave eigenvalue. To the extent that the 2D 
Hubbard model gives an appropriate description of the 
cuprates, this will indicate how well ARPES results can 
be combined with inelastic neutron scattering results to 
provide an estimate of the strength of the pairing inter- 
action in the high-Tc cuprates. A similar analysis using 
Eliashberg equations was applied to the heavy fermion 
superconductor UPtg^i^. 

Fig. [3] shows a comparison of the coupHng U obtained 
from fitting the self-energy (open symbols) and the cou- 
pling obtained from fitting the pairing interaction as de- 
scribed in Sec. IIIII (solid symbols) for various fillings. 
The corresponding d-wave eigenvalue one obtains by us- 
ing the approximative form of the pairing interaction. 




O Zm, <n>=0.95 
-» fit, <n>=0.90 
□■ zm, <n>=0.90 
fit, <n>=0.85 
^ Zfit, <n>=0.85 



O, 

-•- QMC, <n>=0.95 

O RPA, <n>=0.95 

-■- QMC, <n>=0.90 

■a RPA, <n>=0.90 

-♦- QMC, <n>=0.85 

O RPA, <n>=0.85 

02 



FIG. 3: (Color online) The estimate for the coupling strength 
U versus temperature determined from fitting the "exact" 
DCA QMC self-energy with the RPA form [Eq. ([SJ] compared 
to the coupling obtained from fitting the pairing interaction 
with the RPA form [Eq. ([T|] for different site filhngs (n) . 



\U'^X{K - K'), in Eq. (P is shown in Fig. HI Here, as 
a comparison, we also plot the "exact" results for ob- 
tained if the DCA QMC results for the pairing interaction 
T'P'P are used in Eq. ([5]). In Fig.[3]one sees that the single- 
particle estimate of U is smaller than the that obtained 
from fitting the pairing interaction for all temperatures 
and fillings. The result of this is that the corresponding 
eigenvalue obtained in this approximation can be larger 
or smaller than the "exact" eigenvalue, depending upon 
the doping and the temperature. Nevertheless, as shown 
in Table |T1 the errors in the determination of the mean- 
field transition temperature are similar in size to the case 
in which U was determined by a fit that required knowl- 
edge of TPP. 



V. SELF-CONSISTENT DETERMINATION OF 

U{T) AND G{K) 

In the previous sections, we used DCA QMC re- 
sults for the spin susceptibility x{Q) ^^'^ the single- 
particle Green's function G{K), and estimated the cou- 
pling strength U (T) either by fitting the pairing interac- 
tion or the nodal quasiparticle weight. In this section, 
we go one step further and assume a scenario where only 
limited information is available for the single-particle 
Green's function G{K). In this case, we use the approx- 
imate form [Eq. for the self-energy Ti{K) to deter- 
mine the dressed Green's function from the DCA coarse 
grained Dyson equation 



FIG. 4: (Color online) The d-wave eigenvalue versus tem- 
perature obtained from the RPA form [Eq. ^] with U de- 
termined from fitting the self-energy with the form [Eq. ©J 
(open symbols) compared to the eigenvalues obtained from 
the "exact" DCA QMC interaction (solid symbols) for differ- 
ent site fillings (n). 



Eqs. dS]) and dH) are iterated until self-consistency is 
achieved, and the value of tj entering Eq. ^ is again 
fixed by requiring that the quasiparticle weight Z(K,T) 
obtained with this approximation for K = (7r/2,7r/2) is 
the same as that obtained from the "exact" DCA QMC 
self-energy. We find that the estimates one obtains for 
U using this approach are almost identical to the values 
obtained in Sec. lIVK see Fig. [3]) and therefore do not show 
the results. 
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G-i(K,c.„) = [Goi(K + k,c^„) -S(K,c.„ 



.(8) 
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FIG. 5: (Color online) The d-wave eigenvalue Ad versus tem- 
perature obtained from the RPA form [Eq. ([1])] with U and 
G{K) determined self-consistently from fitting the self-energy 
with the approximate form [Eq. (|6}] (open symbols) compared 
to the eigenvalues obtained from the "exact" DCA QMC in- 
teraction (solid symbols) for different site fillings (n). 



Here, Go(fc) is the non- interacting Green's function, i.e. 
Go(k, LOn) = {i^n — Ck)"^ with ek = — 2i(cos -\- cos ky). 



Using these estimates of U{T) and G(fc), one can again 
calculate the d-wave eigenvalue and compare it to the 
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"exact" eigenvalue. The curves with open symbols in 
Fig. [5] show the result for the d-wave eigenvalue obtained 
with this approximation and the curves with solid sym- 
bols display the "exact" result. One sees that the results 
obtained with this approximation are almost identical to 
the results obtained in Sec. IIVI where the "exact" DCA 
QMC Green's function was used. In this case, the tem- 
peratures at which XdiT) = 1 are listed as T^^^ in Table |T] 
and one sees that they are again within 10% to 30% of the 
exact values. Apparently the additional renormalized 
one-loop approximation [Eq. ^] of the self-energy that 
enters the Bethe-Salpether equation [Eq. (O] through the 
propagator G(fc) has only a negligible effect on the d-wave 
eigenvalue. 

We also performed fits of the real frequency spec- 
tra with the corresponding one-loop approximation. On 
small clusters we find it difficult to describe the spectra 
with this simple approximation. On larger clusters, how- 
ever, we have previously shown that the spectra are well 
approximated by the one-loop form for the self-energj*i^. 
Preliminary results on a 16-site cluster indicate that very 
similar estimates of Ubar are obtained as those presented 
in this manuscript. 

VI. CONCLUSIONS 

In conclusion, this work has shown that the pairing 
interaction in the 2D Hubbard model in the parame- 
ter regime appropriate for the cuprates over a range of 
dopings is well described by the spin susceptibility rep- 



resentation, |/7(T)^x(-f^ ~ K')- At low temperatures, 
the coupling strength U decreases with temperature and 
doping. Close to the superconducting transition tem- 
perature, U can be well estimated by assuming that 
the self-energy is determined by the same interaction, 
|C7(T)^x(Q)) S'lid requiring that one has the same nodal 
quasiparticle weight as the "exact" DCA QMC result. In 
practice, one would seek to relate this to the renormal- 
ization of the nodal Fermi velocity measured in ARPES 
studies. Using this approximation to self-consistently 
determine the single-particle propagator that enters the 
Bethe-Salpeter equation in the particle-particle channel 
has only neglibible effects on its d-wave eigenvalue. As a 
result, estimates of the superconducting mean- field tran- 
sition temperature using the spin-susceptibility represen- 
tation of the pairing interaction and the self-energy pro- 
vide a satisfactory way of measuring the strength of the 
d-wave pairing interaction. 
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